Barchan dune's size distribution induced by collisions 
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Large dune fields can overrun for instance infrastructures lying in the way of their motion. Uni- 
directional wind fields and low sand availability give rise to single barchan dunes moving in the 
direction of wind. They can interact by directly exchanging their sand through collisions. This 
kind of interaction play a crucial role in the evolution of barchan dune fields, for instance in the 
selection of a characteristic dune size. Simulations of dune collisions with lateral offset are studied 
systematically and described by general phenomenological rules. Moreover, simulations with a sim- 
plified model without sand flux considering only collisions show that the sizes of sand dunes in such 
a scenario follow a Gaussian distribution with a well defined characteristic size. 



PACS numbers: 

I. INTRODUCTION 

Vast areas of our planet have arid environments where 
non-living matter predominates. The rare vegetation al- 
lows at many places the existence of highly mobile sand 
dunes which move in wind direction up to several decades 
of meters per year. For instance, in Morocco at the lat- 
itude of the Canarian Islands sand is shifted by wind 
from the beaches hundreds of kilometers into the conti- 
nent (see fig. [T]). The unidirectional wind regimes and 
not too high sand coverage lead to the formation of iso- 
lated sand dunes called barchans. They move mostly 
uncorrelated with velocities proportional to their recip- 
rocal height [l[ . The intrinsic instability of barchan dunes 
under an incoming sand flux leads to an increase of the 
largest dunes in the field whereas the smaller ones shrink 
until they disappear @, Hence, the mean size of the 
dunes should grow with the distance from the beginning 
of a field. Nevertheless, in many dune fields the sizes 
saturate. Two mechanisms have been proposed to avoid 
unlimited dune growth: instability of large dunes due 
to changingwind directions Q and collisions between 
dunes [ISMIS- 

Here, we concentrate on the case when a dune collides 
with another one. Due to the size dependence of their 
velocities, the dune that collides onto another one from 
behind must be smaller. This process has been observed 
several times [1, Q but was not understood until recently. 
The large temporal scales of such a process makes it diffi- 
cult to observe the final state after such a collision. Sim- 
ulations were carried out to understand what happens 
if two dunes collide with each other and different results 
were obtained depending on the relative height difference 
between the interacting dunes [7[ ■ The three different fi- 
nal situations of coalescence, solitary wave behavior and 
breeding provide mechanisms to redistribute the dune 
sizes and thus to avoid continuous growth of the dunes 
in a dune field. Here, we present simulations of dune 
collisions with lateral shift in order to extract a general 




FIG. 1: Example of a barchan dune field in Morocco, West 
Sahara. Note that the dune field is divided in corridors along 
the wind direction, where the dune size is roughly uniform. 



rule for the resulting volumes and positions. These new 
rules are implemented into a model for the evolution of a 
simplified dune field where dunes only collide with each 
other. From this model we find that collisions as a mech- 
anism of redistribution of sand capable to select a char- 
acteristic dune size. In fact, dune sizes in such a scenario 
follow a Gaussian distribution. 
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FIG. 2: Snap shots of binary collisions for one volume ratio (r; = 0.2) and different initial dune offsets. For this volume ratio 
no coalescence takes place, and only breeding (a,b,c) and solitary wave behavior (d,e) occurs. 



II. COLLISIONS WITH OFFSET 

This section presents the results of the DUNE model 
that takes into account the calculation of the turbulent 
wind field, the sand flux and the height changes on the 
surface. The DUNE model considers the height profile 
of sand dunes in a unidirectional wind field and calcu- 
lates its impact on the air flow and the sand flux over 
a dune. The dune shape can change and so produce an 
alteration of the atmospheric boundary layer until even- 
tually a final shape forms. The model is explained in 
detail in refs. @, OH HH, EH- Here, we focus exclusively 
on the quantitative results by considering the dunes as 
entities without regarding their particular shapes. 

In the following, the dune sizes will be expressed by 
their width w in order to be able to compare the results 
of the DUNE model in a better way to the results of the 
simulations described in the next sections. The linear 
relations between dune height, width, length and the cu- 
bic one between width and volume are assumed to follow 
the laws observed in measurements [l| and computational 
models [l3j . 

Simulations of the DUNE model are carried out to ob- 
tain the outcome of collisions between two dunes for dif- 
ferent initial lateral offsets 9i and volume ratios r» keep- 
ing the volume of the larger barchan constant at 6000 
m 3 . Open boundary conditions are used with an influx 
equal to the outfiux of the larger dune. We define the 
relative lateral offset between two dunes as 9 = , 
where Y and y are the coordinates of the center of mass 



of the large and the small dune in the lateral direction 
transverse to their movement, respectively. W denotes 
the width of the larger dune. Besides, the volume ratio 
is denned as r = y, where V is the volume of the large 
barchan and v the one of the small one. 

Fig. [2] illustrates different temporal stages of the col- 
lision process for different offsets and one relative size 
r, = 0.2. In general, we distinguish between three fi- 
nal states after a collision: coalescence where only one 
dune remains, solitary wave behavior where the number 
of dunes remains two after the collision (Fig. [5] d,e) and 
breeding where one small dunes leave the horns of the 
larger one (Fig.[Ha,b,c). A similar picture, although less 
rich, was obtained by [Tj| . As can be observed in the fig- 
ures, in the case of breeding only a very small dune leaves 
the horn. This dune generally is not stable. Hence, we 
will neglect the very small dune in breeding and consider 
the output of a collision as consisting only on two dunes, 
a small one that leaves the hybrid state and a large one 
that remains behind. Therefore, breeding will be equiv- 
alent to solitary wave behavior in the following. 

The redistribution of the total volume after a collision 
process is described by the volume ratio r f after the col- 
lision. This volume ratio can approximately be expressed 
by the phenomenological equation, 

77(0*, r*) » 1 - exp (-AQBMn - r o (0i)) 4/3 ) (1) 

valid for n > tq- This condition takes into account that 
there is a minimal relative size robetween the incoming 
dunes below which coalescence occurs [H, Q. This coa- 
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lescence threshold is function of the relative initial dune 
offset and is given by, 

ro(0i) « 0.12 exp(-(^/0.4) 2 ) - 0.05 (2) 

On the other hand, the term A(9i) represents the sen- 
sitivity of the final volume ratio rf to the initial offset 

A(9i) « 10 exp(- (0^0.36) 13 ) (3) 

Fig. [3] illustrates the shape of r/(r*i,#i) compared to 
the measurements (dots) obtained from the simulations 
with DUNE. We analyzed simulations with different off- 
sets 9i for seven different values of r, (Fig. [2] shows the 
simulation results for = 0.2). 

Additionally, the lateral positions of both dunes change 
after the collision. No simple rule could be found and 
thus we do not consider the final offset. 




FIG. 3: Volume ratio 77 after the collision as a function of the 
initial offset 6i and the initial volume ratio Ti due to eq. fl}. 
Dots show the measured values obtained from the simulations. 



dune sizes which depends on the inverse of the square of 
the dune width (Fig. @J , 




where a « 0.45, b « 0.1 and w c are parameter fits to the 
data set. 

The deviation from the scale invariance in the dune 
outflux q out for small dune sizes is expressed by the last 
term of eq. ([5]). However, since the critical width w c is 
of the order of the minimal dune width, the nonlinear 
term is very small and can be neglected. Therefore, we 
consider dune shapes as scale invariant. In this case, the 
mass balance becomes 

«(i->(t-|) w 

where q c = b Q/(l — a) « 0.18 Q denotes the equilibrium 
influx at which the dune volume does not change. How- 
ever, this equilibrium is unstable since there is no mech- 
anism by which barchan dunes can change their outflux 
to adjust it to a given influx. This instability is remains 
after including the correction term for small dunes. 



o 







0.4 




















1 & ° 


0.3 
0.2 
0.1 














D 


0.1 0.2 0.3 0.4 0.5 0.6 













,„ 20 40 60 80 100 120 

w (m) 



III. FLUX BALANCE ON A BARCHAN DUNE 

In order to illustrate the dune size instability let us dis- 
cuss the flux balance equation. The flux balance of a dune 
is given by the difference between the net sand influx Qi n 
and outflux Q ou t- Since both scale with Qw, where Q is 
the undisturbed saturated sand flux on a plane, the mass 
conservation states 

V = Q in -Q out = wQ f^-^) (4) 

where qi n and q ou t are the dune influx and outflux per 
unit length, respectively. 

The simulations of single barchans with a constant in- 
flux using the DUNE model, show that the dune outflux 
scales linearly with the influx with a correction for small 



FIG. 4: Relation between the dune outflux and the dune 
width w. For small widths w < w c the outflux deviates from 
the fit (solid line) which means that the dune has become a 
dome. For higher width however, the outflux relax as 1/w 2 to 
a constant value -solid line. Inset: the dune outflux is a linear 
function of the dune influx -solid line- with a slope smaller 
than unity -dashed line. There are two different regimes, for 
qin < 0.18 Q the outflux is higher than the influx and the 
dune shrinks, while for qi„ > 0.18 Q the influx is larger than 
the outflux and the dune grows. 

In general the flux balance on a dune is not in equi- 
librium. However, assuming that the flux inside the field 
is homogeneous due to diffusion, eq. ^ with an influx 
<lm = Qc leads to the ideal state at which the volumes 
of all dunes in the field remain constant. Under this ap- 
proximation, we consider a simplified model for a dune 
field that only includes binary collisions. 
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IV. THE SIMPLIFIED MODEL 

The system consists of a large number of dunes, char- 
acterized only by their width, which interact exclusively 
through collisions between them. We study the evolu- 
tion of the dune size distribution N(w) in the field in 
order to check if the macroscopic behavior of the system 
approaches an absorbing state. 

We use two kinds of initializations for the system. In 
the first the size of the initial dunes (their width) is ran- 
domly chosen with an uniform distribution between and 
some value w m ax- The second one generates dunes with 
sizes differing only slightly from some previously chosen 
width Wchar (uniform distribution with width two around 
w c har)- The initial number of dunes is 1000 in our simula- 
tions. During a simulation very small dunes are removed 
from the dune field. In order to simulate an entire dune 
field, we consider annealed interactions among the dunes 
which means that a new mutual lateral distance (offset) 
is tossed before each collision. 
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1000 iterations 
— 5000 iterations 
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FIG. 5: Histograms of the size distribution at different times. 
The field is initialized with a distribution strongly peaked 
around w char = 10,20,50, 100 (from left to right). The fig- 
ure shows that the final state is reached very fast (after about 
100 collisions per dune) and a Gaussian size distribution is ob- 
tained. Inset: Half of the points (circles and squares) of the 
figure correspond to the simulations where the initialization 
was with a strongly peaked distribution. The other points 
(diamonds and stars) correspond to an initialization with an 
uniform distribution (w m ax = 10, 20, 50, 100 from left to the 
right). The position of the peak of the size distributions (cir- 
cles and diamonds) as well as their deviation (squares and 
stars) increases with a power law of exponent 1/3. 

For each collision two dunes are randomly chosen from 
the field to collide. This is repeated every iteration as 
many times as there are dunes in the field. If the second 
dune is larger than the first one, no interaction takes 
place. No new dunes are generated but the ones with 
size w < 1 are removed. 

The outcome of a collision between the larger dune 
(upper case variables) and the smaller one (lower case 



variables) is, for simplicity, assumed to depend only on 
the relative mutual lateral distance, the random offset 6i 
and the volume ratio rj = w 3 /W 3 . The new volume ra- 
tio, r/, is cast, according to the measurements presented 
above, through eq. ((T|), and the conservation of the total 
volume of the two dunes, Vtot oc u> 3 + W 3 . As already 
mentioned above, eq. fl} includes only solitary wave be- 
havior. We neglect here the creation of new small dunes 
by breeding. We choose the offset between the two dunes 
to be a random number distributed uniformly in the in- 
terval [0, 0.8] at each collision. The limit of the maximum 
offset of 0.8 is due to the fact that the size of dunes col- 
liding with an offset larger than the distance to the horns 
of the larger one generally is not altered. 

V. RESULTS 

The histograms in fig. [5] show that the dunes form 
Gaussian-like distributions which peak around some 
width w. The results of the model depend on the shape 
of the initial size distribution of the field. The total mass 
of all dunes is conserved with the exception of a negligible 
amount due to small dunes removed from the field. 

The average volume represents the crucial and only 
value to determine the absorbing state of the simulations 
independent from the initial distribution of sizes. The in- 
set in Fig. [5] shows a reciprocal cubic increase of w versus 
the average volume of the field. The distributions also 
become broader for larger w as can be observed in both 
figures. We can calculate the root mean square deviation 
in the corresponding histograms. This deviation follows 
the same simple law as w versus the average volume. 

VI. DISCUSSION AND CONCLUSIONS 

The analysis of collisions of barchan dunes with lateral 
offset showed that some regularity exists in the resulting 
sizes after the collisions. Some simple laws could be for- 
mulated to fit the outcome of such an interaction between 
two dunes. The resulting volumes can be understood to 
depend only on the initial offset and volume ratio. 

We applied these phenomenological rules to a simple 
model which considers the effects of collisions in a dune 
field without considering the sand flux between them. 
The distribution shows a Gaussian-like shape with re- 
spect to the one-dimensional extensions of the dunes. 
The cubic relation between volume and width of a dune 
implies that the volume is not Gaussian-distributed as 
could have been expected. 

We demonstrated that dune collisions could play a cru- 
cial role for stabilizing the size distribution of a barchan 
dune field as already mentioned by ref. [3|. The simple 
model presented here certainly does not resemble real 
dune behavior. Especially the neglection of sand flux 
and the absence of a dependence of the dune interaction 
on spatial dune coordinates hinders us from providing a 
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complete picture of a dune field. The missing mecha- 
nisms need to be incorporated to construct a dune field 
model able to fully understand the effects of collisions 
onto the macroscopic behavior of a dune field. 
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